module gpa

    implicit none
    private

    public :: aaf

    contains

    subroutine aaf(matrix, fa)
    
        real, dimension(3,3), intent(in) :: matrix
        real, intent(out) :: fa
        
    end subroutine aaf

    subroutine gradient(matrix, grad)

        real, dimension(3,3), intent(in) :: matrix
        real, dimension(2,3,3), intent(out) :: grad

        do i=1,3
        do j=1,3
            grad(1,i,j) = 
            grad(2,i,j) =
        enddo
        enddo

end module gpa 
